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Abstract 

We theoretically describe the spin excitation spectrum of a two dimensional electron gas em- 
bedded in a quantum well with localized magnetic impurities. Compared to the previous work, 
we introduce equations that allow to consider the interplay between the Coulomb interaction of 
delocalized electrons and the sd exchange coupling between electrons and magnetic impurities. 
Strong qualitative changes are found : mixed waves propagate below the single particle continuum, 
an anticrossing gap is open at a specific wavevector and the kinetic damping due to the electron 
motion strongly influences the coupling strength between electrons and impurities spins. 



PACS numbers: 75.30.Ds, 73.21.-b, 85.75.-d, 76.50.+g, 76.30.-v 



I. INTRODUCTION 



Collective spin dynamics in dilute magnetic semiconductors (DMS) has recently drawn 
lots of attention.^™ This field provides an insight into the origins of carrier-induced ferromag- 
netism in semiconductors^' and particular features in the spin excitation spectrurrP^ due 
to the presence of two spin sub-systems that are dynamically coupled by Coulomb-exchange 
interaction: that of the itinerant carrier and that of the localized magnetic impurities. The 
transverse spin excitation spectrum has been theoretically found to be composed of three 
types of excitations. These are: two collective spin waves corresponding to itinerant and 
localized spins precessing in phase or out of phase to each other, and single-particle (or 
Stoner-like) excitations of the itinerant carriers.^^. If the DMS is in the ferromagnetic 
state, the in-phase spin wave (IPW) becomes the Goldstone-like mode with an acoustic 
type dispersion responsible for long-range spin order in the ground state. The out of phase 
spin wave (OPW) develops an optical branch, with a zone-center energy determined by 
the strength of Coulomb-exchange interaction between carriers and the spins of magnetic 
impurities. 

Experimental evidence of the entire spectrum in a ferromagnetic DMS like GaMnAs is 
not available. What has been reported so far are features related to the zone-center IPW, 
dominated by the Mn spin precession, its dynamics J 2 * 11 * 12 ! and its ferromagnetic resonance.^ 
We find no experimental data available for the out of phase mode. Indeed, ferromagnetism 
in GaMnAs systems requires a high Mn concentration, which destroys the periodicity of the 
crystal potential and smooths out all optical resonances. 

More insight into the DMS spin excitation spectrum has been gained in CdMnTe doped 
quantum wells (QW), which constitute a clean test-bed system, appropriate to capture 
general properties of the collective spin dynamics in DMS materials. Evidence for carrier- 
induced ferromagnetism has been found in CdMnTe quantum wells doped with holes^. 
When doped with electrons, due to the very low Curie temperature, only the paramagnetic 
phase is available to most experiments. The OPW mode dispersion and single-particle 
excitations have been probed by Raman measurements in the paramagnetic state™ 51 . The 
mixed nature of the IPW and OPW waves has been evidenced in the frequency^ and time 
domairP^l. Neithertheless, there is a lack of a full theoretical description of the spin excitation 
spectrum in CdMnTe QW. Indeed, so far two approaches were followed to describe the spin 
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excitations : in the first one®', the s<i-exchange dynamical coupling between Mn and electron 
spins was considered, but the Coulomb interaction between electrons was dropped out. In the 
second®', the reverse point of view was adopted : spin resolved Coulomb interaction between 
electrons was taken into account, but only the static mean-field sd contribution from the 
Mn was kept to form a highly spin-polarized two dimensional electron gas (SP2DEG). 

This work fills the gap between the two theoretical approaches, by solving the spin dy- 
namical equations in presence of both the sc?-exchange dynamical coupling between Mn and 
electron spins and the Coulomb interaction between electrons. Starting from the full DMS 
Hamiltonian, the approach combines exact commutation rules and standard generalized 
Random Phase Approximation (RPA). We also include the intrinsic damping of the pure 
electron spin waves due to the delocalized character of the electron ^ 19 1 20 1 . We show that the 
introduction of Coulomb electron-electron interaction induces strong qualitative changes in 
the spectrum compared to the approach of RefP^ and that inclusion of the intrinsic damp- 
ing diminishes the strength of the coupling between the two spin subsystems for non-zero 
wavevectors. Generalization of this model to hole systems might be considered : then one 
should take into account the fact that the hole spin states are not isotropic. 

The paper is divided as follows : in Sec. II, we detail the Hamiltonian of the system 
and rewrite it in terms of collective variables, in Sec. Ill, we use transverse spin dynamics 
equations to derive spin response functions, and in the last section, we study the spectrum 
of spin mixed electron-Mn modes. 

II. THE 2DEG DMS HAMILTONIAN UNDER STATIC FIELD 

We consider a QW of width w containing x e ffN unpaired^ Mn spins per unit volume. 
The first subband is populated by electrons per unit surface. The 2DEG-DMS Hamil- 
tonian under the influence of a static magnetic field B = Be z applied in the plane of the 
QW writes : 
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H — Hxin + Hcoulomb + H s _d + Hzeeman (1) 

Hs-d = -a J J J S (r) • M (r) d 3 r 

Hzeeman = g C ^B / / / S (r) ■ BtfV + g M n^B / / / M (r) ■ B<Pr 



where a is the exchange coupling between conduction electrons and Mn spins (a > 0) , 
and g e and gun are normal g-factor of, respectively, conduction electrons and Mn electrons. 
In the convention where /xb > 0, we have g e ~ —1.44 and gMn — 2.00. We have introduced 
two vector operators : S (r) = x 2 (y) ^ §j<5 (r// — IV/J is the 3D electron spin density in a 
splitted coordinates frame r = (r//,y) with r//, the in-plane position and y the out of plane 
coordinate, x (y) is the electron envelope-function of the first subband of the QW. The i 
index accounts for the i-th electron of the 2DEG, its spin | is described by the operator §j 
and its position is v^//. M (r) = . IjS (r — Rj) is the Mn 3D spin density. The j-th |-spin 
Ij of a single Mn impurities is localized on the cation site Hj. In the equilibrium state at 
temperature T, each Mn spin has the average value \IzJ (B, T) , which is the thermodynamic 
average over the five occupied states of the Mn atom d-shell, given by the modified Brillouin 
function^. The 2DEG has the equilibrium spin polarization ( = (n^ — n_\) /n-^o- 

We, now, rewrite the s — d Hamiltonian using the electron (and Mn) spin fluctuations op- 
erators at in-plane wave- vector q, respectively, S q = JjjS(r) e - " 1 ' 1 ~* '< 'd 2 r i jdy = ^ Sie~ iq,I V/ 
and M q = Ije~* q ' Rj '//. Due to the 2D and 3D characters of, respectively, the conduction 
electron and Mn spins subsystems, the electron spin-degrees of freedom naturally couples to 
Mn spin profile weighted by the squared electron wave-function. For later convenience, we 
introduce the following n-profile Mn spin fluctuations operators : 



Mfi =w n X 2n (y) M (r) e-^Hd\ndy = w n J2 ^ (Vi) ^~^ n 



Hence, M q = M q 0) . M q " ) are vector operators verifying the following commuting relations 



i^,X:i?l> ( 2 ) 
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where a, (3, 7 = x, y, z and e a ^ n is the Levi-Cevita tensor. It follows 



H 



a 



s—d 
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^0 Sq • 



(3b) 



In Eq. (3a), we have defined the exchange coupling constant a = a/wL 2 , equilibrium 
averaging () and fluctuation operators : 5 A = A — l^A t 
Finally, we get : 



H s ^ d + H Zeeman = Z e (B) • 5S Z ^ =0 + g M n^ B B ■ SM^ =0 + K ■ 6M^ =0 (4a) 

- a 1 - (£ + , q=0 • M a , q=0 - 5_, q=0 • M£U) (4b) 

- 5*4q=o • <^i,q=o (4c) 

- a £ S q • M« (4d) 

q^O 

with the total bare Zeeman energy of conduction electrons given by : 

Z e (B) = «7i^ Mn (iz) (B, T) - \g c \ fi B B (5) 



where 71 = L y 2 (y) dy is the probability to find the electron in the QW, Nun = 
x e ffNowL 2 is the number of Mn spins available in the QW. In Eq. ([5]), we evidence the 
"Overhauser shift" A and the normal Zeeman contribution, opposite to A (the sd contri- 
bution). Indeed, as g e < 0,(?Mn > and a > 0, the Mn spins are anti-parallel to the field, 
thus, the sd coupling aligns the electron spins anti-parallel to the field, while the normal 
Zeeman aligns the electron spin parallel to the field. When the sd coupling dominates over 
the normal Zeeman contribution and B > 0, both Mn and electron spins are J, '-(^h^ < 
and the 2DEG spin polarization degree ( is also negative. 



The first line, Eq. (4a) gives the mean-field "effective Zeeman" Hamiltonian, where no- 



dynamical coupling between electrons and Mn spins appears. The mean-field Hamiltonian 



naturally introduces the "Knight shift", due to the equilibrium electron spin polarization 
(<S'z,q=o)o = u 2dL 2 (/2, which shifts the Mn spin precession energy: 



K = -a (S , Zjq=0 ) 



1 a 
2w 



™2D |CI 



(6) 



Eq. (4b ) and Eq. (4d ) give rise to a first order sd— dynamical coupling between transverse 



spin degrees of freedom and induce spin-mixed electron-Mn modes of precession. They also 
contain higher orders correlation terms which have been indentified to be responsible for a 
damping 5 . The effects of the correlations contained in Eqs. (4b)-(4d) are out of the scope 
of the present work and will be neglected when they appear. 



III. TRANSVERSE SPIN DYNAMICS 



A. SP2DEG dynamics without the s-d dynamical coupling 



In this paragraph, we take into account only the first line of Eqs.(4a)-(4d). Hence con 



duction electron and Mn dynamics are independent. The 2DEG is polarized by the static 
exchange field of Mn spins, this forms a spin-polarized 2DEG (SP2DEG) as described in 
RefP^. The Hamiltonian which rules the electron dynamics in the SP2DEG reduces to : 



HsP2DEG — H Kin + Hcoulomb + Z e (B) S z cl=0 (7) 

Introducing the the electron creation-anihilation operators, a spin-flip single particle excita- 
tion (SF-SPE) is described by a single electron-hole pair operator c£ ^c^, and electrons 
spin- wave operators, introduced above, are given by S +M = S Xt<i + iS y<cl = Xlk c k-qt Ck >^" 



Let's notice that 



*S + ,q> C k-q',t Ck '^ 



0, such that collective and single particle modes are 



not intrinsically coupled. In the following we will use exact commutation rules to write 
equation of motions for these normal modes of the SP2DEG. We will note ( dA/dt ) 



AH. 



SP2DEG 



J 2DEG 

jih the time derivative of A related to Hsp2Deg only. Further, we will make 



use of linear response theory to derive quantities like : 



A: B 



— lim 



A(t),B 



-itot—et 



dt 



(8) 



which gives the linear response of an observable A to a perturbation coupled linearly to 
B in the considered HamiltoniarP^. 
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A; 13)) has the following equations of motion : 



A: B 



i 

UJ 



MB 



i 

UJ 



MB 



1 

l 



A, B 
A, B 



(9a) 
(9b) 



where A 



A, H 



We will note ((MB 



/iH is the time derivative of A related to the considered Hamiltonian. 

2DEG . ( " 

when A is replaced by ( dA / dt 



2DEG 



1. Single particle modes dynamics 

The kinetic Hamiltonian Hxin = Ylk a ^ c ka c K^ an d the mean- field Zeeman Hamiltonian 
conserve the single-particle modes : 



"k-q,t Ck ^' ^Kin + Z e (B) S Z)<i= o 
(_E k — -Ek-q — Z e ) C k _ q ^Ck,4. 

But the Coulomb Hamiltonian H Cou iomb — \ Yl ^q c k+ q ,a c k'-q,a' c kV'Ck l( T couples a 

q^0,k',o-',k,<T 

single particule mode to multi-pair modes having a spin +1 : 



(10) 



|^ C k— q^ki-l-' Coulomb — 

+ £ V * ( C k-q+q', t C k'^ C k'-q', CT C k,4) 
qV0,k',<r 

~ Yl V *' ( C k-q,t C k',- C k'-q',a C k-q',J 
qV0,k',(T 

where Vq/ is the space Fourier transform of the bare Coulomb interaction. It follows 
from Eq. (10) that Hcouiomb does not conserve the SPE motion, but introduces an infinite 



hierarchy where a single electron-hole pair of the Fermi sea (a SPE) couples to multiple 
pairs having the same global spin. Approximations can be made : in the rhs terms of Eq. 



(10), some conserve the SPE motion and renormalize it, others introduce a scattering effect, 



the so-called spin-Coulomb drag^I, and can be described by an electron-electron scattering 
time^r e _ e . The former consists in making the random phase approximation (RPA) on single 



mode dynamics^, i.e., keeping in Eq. (10) only terms which can be written as a product of 
a SF-SPE with an occupation number n^, and replacing by its average value (n-k >a ) . 
Then the Coulomb factor Vq> has to be replaced by a local field factor G xc which accounts 
for the effective dynamical exchange-field produced by other electrons^ (a part of what 
has been suppressed in making the RPA). Adding a damping rate rj = h/r e - e due to the 
scattering leads to the SF-SPE equation of motion that we will use in the following : 

ih (^ C k-q,t C k,^ ^ G = (#k - #k-q - Z e + 2G XC (4,q=o) o ) 4-^ C *,l i 11 ) 
+ G XC ((n^) Q - (^k-q,t) ) ^+,q - ^ C k-q,t C ia 



Eq. (11) evidences the renormalized Zeeman energy, i.e., the spin-flip energy of single 
electrons : 



Z* — Z e — 2G XC (^S Z) q=o^ (12) 

Compared to the bare Zeeman energy Z e , Z* is enhanced by Coulomb-exchange between 
spin-polarized electrons, a phenomenon linked to the spin-susceptibility enhancement!^. 
Each SF-SPE is characterized by two wavevectors k and q. At q = 0, SF-SPE are de- 
generate to Z*. When q ^ the degeneracy is lifted by the kinetic spread of velocities which 
depend on the initial momentum k. 



2. Collective modes dynamics 

Along the SF-SPE, the above spin polarized SP2DEG develops collective modes, the so- 
called spin-flip waves (SFW). SFW dynamics are described by the S+ )q operators dynamics. 
As Hcouiomb conserves the macroscopic spin, it follows : 



• ,q ; ^-Coulomb 



+ Z e {B)S ZM=0 =-Z e {B)S. 



(13) 



m 



But, the kinetic Hamiltonian couples collective states to the spin current J +q 
7 J2k (k— c k-qt Ck ^ carr i e d by single particle states : 



S-L. n, H 



,qj Kin 



foq ■ J 



(14) 
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Finally the equation of motion of collective SFW modes writes 



^+,q) = iw e £+, q - iq ■ J +iq (15) 

, at J 2D EG 

where u e = Z e (B) /his the frequency of the Larmor's electron mode. 

One is left with evaluating the spin-current dynamics to find the 2DEG electron spin 
waves. The spin current evolution is dominated by single particle states dynamics as Hcouiomb 
does not conserve J + q and destroys the coherence between the single particle objects com- 



posing J+, q . As seen from Eq. (11), the exchange field produced by the spin fluctuation 
SV q drives the spin current. The interplay between the spin-current dynamics and the S +jCL 
dynamics then determines both the SFW dispersion and its damping. The relevant spin- 
current response is the transverse spin-conductivity^ which links the spin current to the 
gradient of the exciting exchange field : 



where () is the expectation value at frequency u. The spin-conductivity has an imaginary 
part originating from the damping of SF-SPE, intrinsically due to Hcouiomb or any source 
of disorder acting on transverse spin degrees of freedom. Consequently, the real part of the 
spin conductivity determines the SFW dispersion, while the imaginary part determines its 
damping. It is worth noting, that the spin wave damping originates from the kinetic motion 
of the conduction electrons and from the topoly of the conduction band, a 2D parabolla. 
In a Luttinger liquicP^, the conduction band is linear in k and ID, thus Eq. (14) conserves 



the macroscopic spin. It breaks the coupling between S+ t q and SF-SPEs, which are coupled 
to charge degrees of freedom by Hcouiomb- This property is at the origin of the well known 



spin-charge separation^ occuring in Luttinger liquids. Injecting Eq. (16) into Eq. (15) and 



solving it in the frequency domain for long wavelength (q <C kp) leads to : 

^£f,q = *W ff 5+, q (17) 
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with u q a complex pulsation: 



Reuiq = co e — q 2 G xc lim Recr + (18a) 
Imuig = q 2 G xc lim Ihkx + (18b) 

q— >Q,uj— »0 

In the following we note a + = lim 9 ^. ,u^o Imcr + , the imaginary spin- conductivity. It 
was calculated in RefP^ and some corrections were added in Ref.^ which gave also an 



experimental evidence of the kinetic damping law found in Eq. (18b). We highlight that 
these q 2 laws are valid in the longwavelength limit when the SFW propagates far from 
the SF-SPE continuum (see RefP^). When close to this continuum, the stronger coupling 
with SF-SPEs introduces corrections to the above laws and one should better replace the 
dispersion law with the pole appearing in the transverse spin susceptibility (see RefP^) which 
will be derived in the next paragraph. 

3. Electron spin- susceptibility 

The transverse spin-susceptibility, defined by the ratio of the expectation value \^S +M 
to the perturbing potentiel g e /XB& +!qaJ , where &+ iqu , is the amplitude at the same pulsation uj 
and wavevector q of the exciting magnetic field, is given by : 

/ / - \ \ 2DEG 

X+ (q,u;) = ((5 + , q ;5_,_ q )> (19) 



Straightforward calculations using the equation of motion Eq. (11) lead to : 



.NHfeUr- .^t) (20) 



18 . 



where we have introduced the transverse Lindhardt-type respons< 

II+ t (q,a;) = y - (21) 

A comparison between the above spin-susceptibility expression and the one given by local 
spin-density approximation^, gives the expression of the local field factor G xc : 

2 1 0E XC 

G xc = g —— (22) 
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where E xc is the exchange-correlation part of the ground state energy^. 
SFW appear as poles of x+ (q^) > one finds in the long wavelength limit, another ex- 
pression for Re u q : 



i z e h 

Reu = uj e — —- — 

q \(\Z*-Z e 2m* 



q 



(23) 



Alternatively, if one uses the approximated equation of motion Eq. (15), one finds the 
spin susceptibility in the long wavelength limit : 



X+ (q,wj 



2(5 



2,q=0 



HbJ — hbJn 



(24) 



B. Transverse spin dynamics equations with s-d dynamical coupling 



Now, we keep lines (4b) to (4d) in the s-d Hamiltonian, and we reconsider collective 
transverse spin dynamics. In the following, the derivative dA/dt = A,H/ih takes into 



account the coherent coupled dynamics due to lines (4b) and mdh in the s-d Hamiltonian, 



(i) 



but reduced to first order terms : higher order correlation terms like J^ q , 5S z ^ +q i ■ M_j_ 
have been dropped. 



1. Electron dynamics 



We find : 

where 5 +iq = iu q S +tCl . Compared to the SP2DEG dynamics, the s<i-dynamical coupling 



= s+ « ~ n KM ^ (25) 



adds the second term of Eq. (25) which is a coherent coupling with Mn transverse degrees 



of freedom. One key feature is that the collective electron motion naturally couples with 
Mj 1 ^ Mn-modes, a Mn precession having a profile, out of the QW plane, following the 
electron probability distribution. We are left with deriving the equation of motion for these 
Mn-modes. 
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2. Manganese dynamics 



We obtain the first order equation of motion for Mn spins : 



d _ i n ,. Dl/W i i 1^ »>("+!) z 



= -g Mnf M B BM™ + .rKM^ - ^A„, +1 S + , q (26) 



rft + ' q r JW n a 

where we have introduced n-profile Overhauser shifts : A n = a 



1 ' 



7n/7lA 



with 7„ = w ra 1 J" x 2n (y) dy. The important features are the second and third terms in 



Eq. (26). The later couples the Mn-precession with collective electron modes. The former 
couples a n-profile Mn mode to a (n + l)-profile mode, because this coupling is mediated 
by the 2DEG. Thus, the Mn-dynamics is given by an infinite serie of equations. This is 
a consequence of the 3D nature of the Mn dynamics. A variable like M+™ q describes an 
oscillation propagating in the plane with a rigid profile in the normal direction, but the out 
of plane degree of freedom is restored by the possibility for Mn spins to build modes which 
are combinations of M+ n q resulting in different out of plane profiles^. Obviously, the M+l 
are not independent variables because they don't correspond to orthogonal out of plane 

~ (n) 

profiles. Solving the serie of infinite equations requires a projection of M\_ q over a set of 
modes with orthogonal profiles as it was carried out in Ref.^. Along with modes having a 
strong mixed nature (electron-Mn modes), we then expect to find a high number of modes 
having essentially a Mn character, but with orthogonal profiles (Mn modes). The number 
of Mn modes has to be consistent with the initial number of degrees of freedom present in 
the system. Ref.^ found a high number of Mn modes branches which were separated by 
energies of the order of 0.1 A. However, in the experimental data of RefP, only one branch 
of these Mn modes was apparent. It appears then, that the set of modes chosen in RefJ-^is 



not the most appropriate to describe properly all the modes contained in Eqs.(25)-(26), at 
least in the vicinity of the anticrossing gap (see below). Anyway, this point requires further 
developments out of the scope of the present study. Indeed, we are particularly interested 
in discussing mixed electron-Mn modes which are strongly coupled to electrons rather than 
modes specific to the 3D nature of the Mn dynamics. We can remark that the coupling 
between M+ n q and M^ 1 ^ has a strength given by K, which is very small compared to 



A„ due to the ratio 



^> 1. Hence, considering only modes strongly 
coupled with electron modes is reasonable. As electron modes are naturally coupled to M+L 
modes, we will consider the dynamics for these ones only by cutting the infinite serie with 
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an nomothetic approximation 



^ = (72/71)^ (27) 
Consequently the set of coupled electron-Mn equations reduces to : 



it 5 + , q = z^5 + , q -^™+,q (28a) 



-Mfi = ^Mti^q - ^ A 25 + , q (28b) 



where 



(gun^BB + K^J /h (29) 



is the natural precession pulsation of the free M^ q mode. 



IV. MIXED MN-ELECTRON SPIN WAVES 



A. Spin susceptibilities 



To find the dynamically coupled modes, we will derive the electron spin susceptibility 



with help of equations of motion (9a) and (28a)- (28b) 



From Eqs.(9a) and (28a), we first get : 



iluJ \' S °/ 



then from Eq. (28b), we get 



A2 

tVuJ 



hence, 



-A 5 



5'+,qi 3-,-q. 



hid — fojJMr 



which finally leads to : 
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1 1 u fiw- hw Mn - 9L J l x+ (q, w) 



and 

M«;M« q = - \ (31) 

/ / w feu - to Mn — (q, w) 

Consequently, e-Mn mixed spin excitations appear as poles of the above responses, i.e., 
are zeros of the propagator : 

huj - hujun - —^x+ (q, w) (32) 

with x + (q,w) being the spin-susceptibility of the uncoupled SP2DEG described in Section 

mm 

We can understand the above equation as follows. Consider the Mn point of view ; 
in the presence of the SP2DEG, the precession frequency of Mn spins is shifted from the 
normal precession [g Mn \x B B jYi) by two quantities : a blue shift due to static exchange field 
with spin polarized electrons (K/h) and an additional shift due to the dynamic change of 
the electron spin-polarization. The later is induced by the Mn precession itself. Finally 



Eq. ( 32 ) describes a recursive closed loop where : Mn transverse precession induces electron 
transverse precession proportional to A2X+ (q, 1 ^) , this dynamically changes the electron spin 
polarization which in turn shifts the Mn precession frequency by an amount &A2X+ (qA>) • 



Finally, in dropping correlation terms given by Eqs. (4c) and (4d), one finds a collective 
behavior where electrons and Mn respond adiabatically to the dynamical perturbation from 
the opposite spin-subsystem. 

Similar expressions for the coupled modes propagator have been obtained in previous 
works. To our knowledge it was first derived in RefPfor bulk DMS, and more recently using 
a spin-path integral approach in DMS quantum wells with electrons^^ or bulk DMS with 
holes 8 . However, none of these works did include the influence of the Coulomb interaction 



between carriers. Instead of Eq. (32), they resulted in the following propagator : 



frhJ - fthJMn ^~ n ^t (q> W ) ( 33 ) 

where the electron spin-susceptibility was replaced by the non-interacting single-particle 



response introduced in Eq. (21). Introduction of Coulomb interaction results in strong 



14 



qualitative changes in the spectrum which have been partially adressed in Ref.^ and will be 
detailed below. 

B. Homogeneous modes 



Mn-like 




4 6 
Magnetic field B (T) 
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FIG. 1: (Color online). Zone center electron-Mn modes. Dotted lines are the uncoupled electron 



(upper curve) and Mn (lower curve) modes. Full lines are the solutions w^L °^ Eq.(34). Out of 



the resonant field Br where the modes anticross, branches have electron or Mn character. Sample 
parameters are, x e ff = 0.23%, T = 2 K, w = 15CL4 and = 3.1 x 10 cm~ 2 . 



Eq. (33) was used to successfully fit the experiment of Teran et al. (Ref.^J) where ho- 



mogenous modes (q=0) were probed and shown to experience an anticrossing at a magnetic 
field Br such that UMn — w e . It is a consequence of the Larmor's theorem that the ho- 



mogeneous electron mode behaves as if electron were not interacting (Eq. (15) for q = 0). 



Indeed, setting q = in Eq. (32), leads to the homogenous precession modes equation : 



(34) 



which solutions are real : 
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Figure [T] shows the magnetic field dependence of these modes and the gap opening at the 
resonant field 5r. The upper branch ojg =0 (resp u~ =0 ) has an electron character (resp. Mn 
character) when B < B-& and vice- versa for B > B^. The amplitude of the homogeneous 
anticrossing gap Aw g=0 = y/4:KA 2 denotes the strength of the dynamical coupling between 
the two spin subsystems. Detailed discussions on this gap have been given in Ref.™ and 
RefP. In particular, RefP identifies the anticrossing Mn mode as M+ \ consistent with the 



homothetic approximation of Eq. (27) used here. The anticrossing gap was found to be 

Acj g=0 = \j AKA 2 — ^zf~) " wnere if" * s the damping rate of the homogeneous uncoupled 
electron mode, a quantity that we have neglected here in dropping electron-Mn correlation 



terms contained in Eqs.(4b)-(4d). In RefP, was estimated from measurements of the 
electronic spin wave damping at q — 0. A rigorous simultaneous determination of A2 and T 2e , 
lead to the extraction of K from the anticrossing gap and furthermore to the spin-polarization 
degree ( of the 2DEG. Data showed that the so-extracted ( was slightly exceeding the 
prediction^ made for ( in contradiction with other determinations of ( performed in the 
same type of samples^, which showed that the model used to predict (was reliable. A more 
accurate description of the anticrossing gap taking into account the infinite set of coupled 
n-profile mode M+ q equation of motion might overcome this discrepancy. One should also 
mention, that the infinite serie of equations must be cut in order to conserve the initial 
number of degrees of freedom (number of available spins in the system). But finding the 



right number of Eq. (26)-like to be taken depends on how the total number of degrees of 
freedom separates into a number of (quasi-) individual modes^ and a number of collective 
electron-Mn modes. Determining this separation is also an important and interesting issue. 

C. Spin waves 



For q > 0, Eqs(32) and (33) give very different qualitative results as illustrated on Fig. |2| 
Without Coulomb interaction between electrons, uncoupled modes of the electrons are the 
SF-SPE which are degenerate to Z e at q = 0. The sd dynamical coupling introduces two 
additional collective mode : the OPW propagating above the SF-SPE domain with a positive 
dispersion and the IPW propagating below with a negative dispersion. Introducing the 
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Wave vector q ((am -1 ) 



FIG. 2: (Color online). Illustration of the changes introduced by the Coulomb interaction on 
the mixed modes dispersions for the same sample parameters as in Figjlj (a) Solutions without 



Coulomb given by Eq.(33 ) for B = 5.8 T. At q = SF-SPE are degenerate to Z e . The sd dynamical 
shift introduces two propagating waves : above (resp. below) the SF-SPE domain, the OPW (resp. 
IPW) propagates with a positive (resp. negative) dispersion. The dashed line is the uncoupled Mn 
mode (degenerate to HuMn)- (b) and (c) The Coulomb interaction between electrons is included 



as in Eq.(32) of this work. Dashed lines are the uncoupled electron and Mn modes. At q = the 
SF-SPE energy is shifted to Z* . An anticrossing gap opens at the wavevector q# (B) if B < Br. 
In ususal conditions (see RefPj) both the OPW and IPW propagate below the SF-SPE continuum 
with negative dispersions. Dipersions were calculated after setting to zero the kinetic damping 
rate. 



Coulomb interaction between electrons shifts the SF-SPE to higher energies (Z*) and givs 
rise to the collective wave SFW propagating below the SF-SPE continuum. The SFW is 
further coupled to Mn modes through the sd interaction. An evaluation of the coupling 
between SF-SPE and Mn modes was given in Ref.^ and found to be negligible. Thus the 
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Coulomb interaction introduces a shift between the SF-SPE and the SFW energies, and the 
later is further shifted by the sd dynamical coupling. In realistic conditions, it was shown in 
RefP^that the Coulomb shifts dominates over the sd dynamical shift. Hence, when Coulomb 



interaction is taken into account, Eq. (32), except under unrealistic conditions, gives rise to 
two spin wave modes propagating below the SF-SPE continuum, the IPW and OPW. An 
anticrossing gap opens at a specific wavevector q-& (B) given by : 



2tt? * 

q R (B) = J |C| (Z*/Z e - 1) — (u e - u Mn ) (36) 

Note that q R (-Br) = and that u) 9R = ujun + ^rG xc o" + . If q R (B) > 0, compared to the 
homogenous gap, the anticrossing gap at q-& is dramatically reduced by the kinetic damping 
of the electron wave and is given by : 



< - *>* = V 4KA ^ 2 - < ( 37 ) 

where r) qR = q^G xc a + . We note that the kinetic damping is the only one considered here. 



Other sources of damping, as e.g., the ones dropped in Eq. (4d), will of course further reduce 
the amplitude of the gap. 

Fig. [3] illustrates the variation of qn (B) with the magnetic field. It is always smaller 
than q m = kpi — kp^, the wavevector delimiting the window where the SFW propagates^. 
Overlaid in Fig. |3j are the anticrossing gap at qn and the corresponding damping rate 
?7g R . In the absence of the kinetic damping, the gap would be given by \/4KA 2 . One sees 
the dramatic effect of this intrinsic kinetic damping, which kills the gap outside a very 
narrow range of magnetic fields. As coupling between spin waves of the electron and the 
Mn spin systems is responsible for the appearance of the carrier induced ferromagnetismP, 
we might conclude that the above disappearance of the gap diminishes the possibilities for 
ferromagnetic transitions with complex order (out of q = 0). 

The disappearance of the gap is illustrated in Fig. |4]by comparing the dispersions obtained 



from the zeros of the propagator in Eq. (32) in the presence or absence of the kinetic 



damping. In the presence of the damping, the solutions have a non-zero imaginary part 
for q > 0. The corresponding damping rate is plotted in the lower insets of Fig. |4} It is 
well known that when the frequencies of two coupled oscillators anticross each other, their 
corresponding damping rates cross themselves. Clearly for B < 5.7 T, the mixed modes do 
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FIG. 3: (Color online). Left axis : variation of the wavevector gp, (B) and the wavevector q m = 
kp^ — kp^. Right axis : variation of the anticrossing gap and the kinetic damping rate ij qR = Imu) qR . 
The anticrossing gap is killed by the damping rate when the later is of the same order of magnitude 
as y/AK&2- A break in the horizontal axis scale has been introduced to zoom the region close to 
B ^ Br. To calculate the damping rate, we have used a typical SF-SPE scattering time r = 2ps 
instead of r e _ e (~ 150ps) to match the experimental conditions of RefP^l 

not anticross at any q and each branch conserves its former character, Mn-like or electron- 
like. On the contrary, for 5.7 T< B < Br, the modes anticross at q-R, and the OPW 
transfers the kinetic damping (q 2 law) of the SFW to the IPW when q > q&. It is worth to 
note that this q 2 law for the IPW damping rate was also found in GaMnAs compounds in 
the ferromagnetic stated 
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FIG. 4: (Color online). Illustration of the progressive disappearance of the anticrossing gap when 
B goes away from Br. Continuous lines (open symbols) are dispersions calculated in absence 
(presence) of the kinetic damping. Dashed lines are the uncoupled modes. Upper insets : zoom on 
the anticrossing gap region. Lower insets : variation with q of the damping rate of both the OPW 
(dashed line) and IP W (straight line). 



V. CONCLUSION 



In conclusion, we have introduced equations for the spin dynamics in a test-bed diluted 
magnetic system that allow to take into account the interplay between the Coulomb interac- 
tion dynamics and the sd dynamical coupling between the electrons and the localized spins. 
We have shown how the Coulomb interaction introduces strong qualitative changes : the 
mixed electron-Mn modes propagate below the SF-SPE continuum and an anticrossing gap 
is open for a given range of magnetic field. Because of Coulomb interaction, the intrinsic ki- 
netic damping due to the electron motion is always present (the SF-SPE scattering time can 
not be longer than r e _ e ), this damping kills the anticrossing gap outside a very narrow range 
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of magnetic fields. Our calculations illustrate also how this kinetic damping is transfered to 
the IPW, a phenomenon found in GaMnAs compounds. 
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